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Abstract 

We have developed a multiple image method to compute the interparticle 
force for a polydisperse electrorheological (ER) fluid. We apply the formalism 
to a pair of dielectric spheres of different dielectric constants and calculate the 
force as a function of the separation. The results show that the point-dipole 
(PD) approximation errs considerably because many-body and multipolar in- 
teractions are ignored. The PD approximation becomes even worse when the 
dielectric contrast between the particles and the host medium is large. From 
the results, we show that the dipole-induced-dipole (DID) model yields very 
good agreements with the multiple image results for a wide range of dielectric 
contrasts and polydispersity. The DID model accounts for multipolar inter- 
action partially and is simple to use in computer simulation of polydisperse 
ER fluids. 

PACS Number(s): 83.80.Gv, 82.70.Dd, 41.20.-q, 02.60.Nm 
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I. INTRODUCTION 



Polydisperse electrorheological (ER) fluids have attracted considerable interest recently 
because the size distribution and dielectric properties of the suspending particles can have 
significant impact on the ER response Real ER fluids must be polydisperse in nature: 
the suspending particles can have various sizes or different permittivities. In a recent paper 
0, we investigated the case when the suspending particles are of different sizes. In this 
work, we extend the study to that of different dielectric constants. 

The point-dipole approximation [^] is routinely adopted in computer simulation [^,|^ 
because it is simple and easy to use. Since many-body and multipolar interactions between 
particles have been neglected, the predicted strength of ER effects is of an order lower than 
the experimental results. Hence, substantial effort has been made to sort out more accurate 
models Recently, we have developed a multiple image method and an integral equation 

approach to compute the interparticle force. In particular, we proposed a dipole-induced- 
dipole (DID) model for efficient computer simulation of polydisperse ER fluids p[. 

Poladian P] claimed that the multiple image method can be used to calculate the dipole 
moment of a pair identical dielectric spheres in an applied electric field. In Ref. 0, we 
generalized the multiple image method to a pair of dielectric spheres of different sizes. We 
showed that the generalization yields a reasonable approximation when the spheres have 
a large dielectric constant. The multiple image method was widely adopted [llO|-|T^. The 
approximation is reasonable because in ER fluids, the dielectric constant of the particles 
can be much larger than that of the host fluid. However, the results for low contrast are 
questionable. 



In fact there is a more complicated image method for a dielectric sphere |]13[ , which gives 
the exact image dipole moment of a dielectric sphere that placed in front of a point dipole. 
We thus modify the multiple image formula. The results of the improved formula agree 
with the numerical solution of an integral equation method ||T^JT5[| even when the dielectric 
contrast of the spheres is low. 
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In this work, we extend the multiple image method to compute the interparticle forces 
for a polydisperse mixture of dielectric spheres of different dielectric constants. The DID 
model will be compared with the Klingenberg's empirical force expressions [0. 



II. IMPROVED MULTIPLE IMAGE METHOD 



Here we briefly review the method and extend the method slightly to handle different 
dielectric constants. Consider a pair of dielectric spheres, of radii a and 6, dielectric constants 
ei and e'^ respectively, separated by a distance r. The spheres are embedded in a host medium 
of dielectric constant t2- Upon the application of an electric field Eq, the induced dipole 
moment inside the spheres are respectively given by: 



(1) 



where the dipolar factors /3, j3' are given by: 

(3 



ei - €2 



ei + 2e2 ' e[ + 2e2 

From the multiple image method p|, the total dipole moment inside sphere a is: 
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where the subscripts T (L) denote a transverse (longitudinal) field, i.e., the applied field is 
perpendicular (parallel) to the line joining the centers of the spheres. Similar expressions 
for the total dipole moment inside sphere h can be obtained by interchanging a and 6, as 
well as (3 and (3' . The parameter a satisfies: 

r"^ — — b"^ 
cosh a = ■ 



2ah 

In Ref. 0, we checked the validity of these expressions by comparing with the integral 
equation method. We showed that these expression are valid at high contrast. Our improved 
expressions will be shown to be good at low contrast as well (see below). 



The force between the spheres is given by [|T6 

Ft = ^^{PaT+PbT), Fl = ^^{PaL+PbL)- (5) 

For mono disperse ER fluids, Khngenberg defined an empirical force expression |^: 

^ = (2F|| cos^ e-F^ sin^ ^)r + Fr sin 260, (6) 



FpD 

being normalized to the point-dipole force Ppp, = —Sp'^/r'^, where F\\, F± and (all tending 
to unity at large separation) are three force functions being determined from numerical 
solution of Laplace's equation. The Klingenberg's force functions can be shown to relate to 
our multiple image moments as follow (here a = b, (3 = (3' and Pa = Pb)'- 

where pL = Pl/FpdEq and pp = Pt/FpdEq are the reduced multiple image moments. We 
computed the numerical values of these force functions separately by the approximant of 
Table I of the second reference of Ref . and by Eq. (^ . 

In Fig.l, we plot the multiple image results and the Klingenberg's empirical expressions. 
We show results for the perfectly conducting limit {(3 = 1) only. For convenience, we 
define the reduced separation a = r/{a + b). For reduced separation a > 1.1, simple 
analytic expressions were adopted by Klingenberg. As evident from Fig.l, the agreement 
with the multiple image results is impressive at large reduced separation a > 1.5, for all 
three empirical force functions. However, significant deviations occur for cr < 1.5, especially 
for For a < 1.1, alternative empirical expressions were adopted by Klingenberg. For 
F±, the agreement is impressive, although there are deviations for the other two functions. 
From the comparison, we would say that reasonable agreements have been obtained. Thus, 
we are confident that the multiple image expressions give reliable results. 



III. DIPOLE-INDUCED-DIPOLE MODEL 



The analytic multiple image results can be used to compare among the various models 
according to how many terms are retained in the multiple image expressions: (a) point- 
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dipole (PD) model: n = 1 term only, (b) dipole-induced-dipole (DID) model: n = 1 to 
n = 2 terms only, and (c) mult ipole- induced- dipole (MID) model: n = 1 to n = oo terms. 

The multiple image expressions [Eqs.(3)-(6)] allows us to calculate the correction factor 
defined as the ratio between the DID and PD forces: 

Jf^ ~ ^ " (r2 - 62)4 ~ (^2 _ q2)4 + (r2 - a2 - 62)4 ' U 

~ ^ (r2 - 62)4 + (^2 _ ^2)4 + (r2 - a2 - 62)4 ' 

where = SpaoPbo/^'^ and = —QpaoPbo/^'^ are the point-dipole forces for the trans- 
verse and longitudinal cases respectively. These correction factors can be readily calculated 
in computer simulation of polydisperse ER fluids. The results show that the DID force de- 
viates signiflcantly from the PD force at high contrast when /3 and (3' approach unity. The 
dipole induced interaction will generally decrease (increase) the magnitude of the transverse 
(longitudinal) interparticle force with respect to the point-dipole limit. 

In a previous work , we examine the case of different size but equal dielectric constant 
(/3 = (3') only. Here we focus on the case a = b and study the effect of different dielectric 
constants. In Fig.2, we plot the interparticle force in the transverse fleld case against the 
reduced separation a between the spheres for f3 = 1/3 and various (3'/f3 ratios. At low 
contrast, the DID model almost coincides with the MID results. In contrast, the PD model 
exhibits signiflcant deviations. Similar conclusion can be drawn from the longitudinal fleld 
case as in Fig. 3. It is evident that the DID model generally gives better results than PD for 
all polydispersity. 

At higher contrast, the DID model still agrees with the MID model, except at close 
encounter. In Figs. 4 and 5, we plot the force in the transverse and longitudinal fleld cases 
against the reduced separation a. It is evident that the DID model generally gives better 
results than PD for all polydispersity. For the longitudinal fleld case, the DID model agrees 
with the MID model for a > 1.2, except at close encounter where the MID force diverges as 
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CONCLUSION 



In summary, we have used the multiple image to compute the interparticle force for 
a polydisperse electrorheological fluid. We apply the formalism to a pair of spheres of 
different dielectric constants and calculate the force as a function of the separation. The 
results show that the point-dipole approximation is oversimplified. It errs considerably 
because many-body and multipolar interactions are ignored. The dipole-induced-dipole 
model accounts for multipolar interactions partially and yields overall satisfactory results in 
computer simulation of polydisperse ER fluids while it is easy to use. 
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Figure Captions 

Fig.l: Comparison of the multiple image results with Klingenberg's empirical force expres- 
sion. 

Fig. 2: Interparticle force for transverse field, /3=l/3 while (3' / (3 ranges from 1.0 to 1.2. 

Fig. 3: Same as Fig. 2. Force for longitudinal field. 

Fig.4: Force for transverse field, /3=9/ll while (3' / (3 ranges from 1.0 to 1.2. 

Fig. 5: Same as Fig.4. Force for longitudinal field. 
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Fig.2: Force for Transverse Field, (5—1/2, 

P'=P P'=1.05P 





1 , 1 , 1 


- \ 


IVllU 


^ \ 


X X DID 




PD 




1,1,1 



1.0 1.1 1.2 1.3 1.4 1.5 
P'=1.1p 





1.0 1.1 1.2 1.3 1.4 1.5 

p'=1.2p 




1.0 1.1 1.2 1.3 1.4 1.5 



1.0 1.1 1.2 1.3 1.4 1.5 



9 



0.00 
-0.01 \- 
-0.02 
-0.03 



-0.05 



0.00 
-0.01 \- 
-0.02 
-0.03 
-0.04 h 
-0.05 



Fig. 3: Force for Longitudinal Field, /9=l/3 
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